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Abstract 

We consider generalized chameleon models where the conformal coupling between 

matter and gravitational geometries is not only a function of the chameleon field (p, but 

bX) also of its derivatives via higher order co-ordinate invariants (such as df^cpd^cp, 0(pj ...). 

Specifically we consider the first such non-trivial conformal factor A{(j)^d^(j)d^(j)). The 

^7r associated phenomenology is investigated and we show that such theories have a new 

^\ generic mass-altering mechanism, potentially assisting the generation of a sufficiently 

^l!^ large chameleon mass in dense environments. The most general effective potential is 

\^ derived for such derivative chameleon setups and explicit examples are given. Interest- 

(y-^ ingly this points us to the existence of a purely derivative chameleon protected by a 

^^ shift symmetry for (f) ^ (p + c. We also discuss potential ghost-like instabilities associ- 

^-_l ated with mass-lifting mechanisms and find another, mass-lowering and instability-free, 

►;>. branch of solutions. This suggests that, barring fine-tuning, stable derivative models 

*k> are in fact typically anti-chameleons that suppress the field's mass in dense environ- 

\—^ ments. Furthermore we investigate modifications to the thin-shell regime and prove a 



no-go theorem for chameleon effects in non-conformal geometries of the disformal type. 
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1 Introduction 

Modified gravity theories (see e.g. [H [2] for reviews) have enjoyed continued interest over 
the past few decades and recent developments have led to a much better understanding of 
scalar-tensor theories in particular. The notion of screening mechanisms - how a light scalar 
degree of freedom (j) can act as dark energy on cosmological scales while being shielded 
in dense environments such as on earth - has turned out to be especially useful in this 
context. Implementations of such a mechanism include the following: 1) The chameleon 
model [31 [1], where a density dependent mass is generated and the field becomes too 
massive for detection in dense environments. 2) Vainshtein screened setups [5l|6l [7j such as 
DGP |8] and Galileon [9]/Horndeski [10] models, where non- linear interactions of lead to 
strongly coupled dynamics. A density-dependent (classical) renormalization of the kinetic 
energy there results in an effectively decoupled scalar in dense environments. 3) Symmetron 



models [HI [12], where a scalar is coupled to matter with a coupling strength proportional 
to the vacuum expectation value of (j). This in turn depends on the ambient density, so 
that the scalar effectively decouples in high-density regions. All these mechanisms reconcile 
the existence of a light cosmological scalar with tight fifth force constraints on solar-system 
scales [13]. 

In this paper we wish to focus on chameleon models and potential extensions thereof. 
Chameleon phenomenology has already been studied extensively [HI [TSl HSl HH UHl UHl 
[2ni[2Il[221l23ll2a[2Sl[2Sl[2Z] for variants of the model introduced by [21 S]. Typically 
such theories are built by universally coupling matter to a metric conformally related to 
the Einstein metric g^^. Making this bimetric structure explicit we here investigate whether 
there are any interesting and qualitatively new implementations of the chameleon mechanism 
that arise when going beyond the simple conformal relationships considered so far. Can the 
chameleon effect occur in new guises for generic 4D effective field theories? In other words, 
we will be asking two questions: What forms can the conformal relation generically take? 
And does it have to be conformal or are there more general bimetric structures that can 
produce chameleonic phenomenology? 

A systematic way of undertaking this investigation is to construct a generic matter metric 
as a function of the Einstein metric g^y, and derivatives of the field c?"'^. By showing how 
what we dub derivative chameleons arise in this framework, this paper aims to illustrate how 
such an approach unveils qualitatively new constructions. More specifically we show that 
derivative chameleons naturally give rise to a new mass-altering mechanism, which changes 
the mass of oscillations around an effective potential minimum. The mass-lifting branch 
may help in ensuring can escape detection in fifth force experiments and consequently 
may be of use in alleviating fine-tuning constraints for chameleon models. However, we 
also show that care needs to be taken in order to avoid ghost-like instabilities for mass- 
lifting solutions. The mass-lowering branch of solutions, on the other hand, is generically 
stable. Importantly the new mechanism works for purely derivative conformal factors too, 
opening up the exciting possibility of having a chameleon mechanism which comes endowed 
with a shift-symmetry in the field — ?■ + c, offering better protection from quantum- 
corrections. Furthermore we discuss modifications to the radial solutions around spherical 
matter sources and modifications to the thin-shell mechanism. We also point out why non- 
conformal geometries are not expected to display chameleon-screening, establishing a no-go 
theorem for so-called disformal geometries. 

The plan for the paper is as follows. In section [2] we review the standard chameleon 
picture and how this gives rise to an effective potential with a large mass for oscillations 
around the effective minimum. Section [3] then introduces the bimetric framework, explicitly 
formulating metric relations which are used to construct chameleon models thereafter. This 
construction takes place in [4], while focusing on conformally related metrics, pointing out 
the phenomenology of our generalized chameleon model, in particular the new mass-altering 
mechanism. In [5] this is followed up by presenting a few simple concrete examples that im- 
plement this new-found mechanism and show how it can be realized for purely derivative 



conformal factors. Modifications to thin-shell screening are then discussed while investigat- 
ing radial solutions around a massive source in section |6] In [7] we go beyond conformal 
relationships, arguing that their natural extension - disformally related geometries - cannot 
produce chameleon phenomenology (they do, however, naturally generate Vainshtein-type 
screening solutions). Finally we conclude inlsl. 



2 Conformal chameleons I: The minimal theory 

Chameleon models |21 E] are typically constructed with the use of two conformally related 
metrics, where the conformal factor is a function of the chameleon field only. In particular 
we may consider an action of the following form 

S = j d'x^ (^R + X- V{<P)] + Sm {A'{<P)g,., ^,) , (2.1) 

where X = — |V^0V^0 is the usual canonical kinetic term for 0, V{(f)) is some arbitrary 
potential and Sm describes the matter part of the action with all matter species \l/j universally 
coupled to the matter metric g^i, = A^{(j))g^y. This universal coupling ensures the validity 
of the weak equivalence principle by design. Also note that the signature of g^i, used here 

is ( — h ++), and hence 2X = 0^ — ( V0) . As usual a matter stress-energy tensor can be 
defined with respect to the matter ("Jordan frame") metric g^^, so that 



f^^ = ^|^. (2.2) 

Since it minimally couples to ^^jy, it is covariantly conserved with respect to that metric 
Vfj^T^'^ = 0. When mapping T^^ into the Einstein frame one finds 

f^u ^ A-%(f))T^'\ (2.3) 

This means T^'^ is not covariantly conserved in the Einstein frame VuT'^^ 7^ 0. Turning our 



attention back to the original action (2.1), one can now work out the equation of motion for 
the scalar obtaining 

n0 = y^ - A^{(f))A,^f, (2.4) 



where the matter stress-energy tensor is defined as in (2.2) and has been contracted with the 
matter metric T = T^^g^^. Specializing to the case of a pressureless, non-relativistic source, 
the only non-vanishing component of the stress-energy tensor Tj^ is Tg = —p. Following 
from (2.3) the energy density of matter in the Einstein frame, p, is given by p = A'^p. As 



a direct consequence of V^T^*^ = we then also find a conserved quantity in the Einstein 
frame p = A^p = A~^p. In terms of this conserved quantity (2.4) can therefore be written 
as 

a^ = V^ + A^^p, (2.5) 
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and one can integrate up to obtain an effective potential for 

VM) = V{<P)+pA{<P), (2.6) 

since p is conserved and independent of in tlie Einstein frame. 

This setup has interesting phenomenological consequences. Suppose we start with a 
runaway potential V{(j)), e.g. V{(j)) = Mp'^'^/cp"' as desirable from the point of view of 
quintessence models [28]. This ensures that, in the limit when we can ignore the matter 
action Sm (a low-density environment with p — )■ in the language set out above), we recover 
a quintessence-like solution to the cosmological constant problem. A slow-rolling light scalar 
field then drives accelerated expansion of space-time. But in regions of high density this 
behavior changes in the following way. Naively does not possess a mass term at all, since 



V{(f)) has no minimum (except for the trivial one at — )■ oo). However, from (2.5) it becomes 



clear that, for non-zero p, a suitably chosen A(0) can result in a Vgff which does have a 
minimum 0min, such that 

Ve//,</,(</>min) = V4>{4>min) + ^,0(</>mm)P = 0, (2.7) 

where the mass of the field m for small oscillations around the minimum 0mm is given by 

"^^ - K;//,9i</.(0min) = K</'</'(0min) + M,0?i(0min)- (2.8) 

This is the essence of the chameleon mechanism: An environmentally-dependent way of 
generating a large mass for an otherwise very light scalar 0. This reconciles a model such 



as (2.1) with fifth force constraints, since becomes too heavy for detection in laboratory 
experiments on earth, yet can act as dark energy on large scales. Figure [T] illustrates the 
chameleon mechanism for a conformal factor of the form v4(0) ~ e''^'^. 

Two final comments are in order here. Firstly our argument has been purely phenomeno- 
logical in that a suitably chosen conformal factor can give rise to the chameleon mechanism 
described, but no arguments for why such an A(0) is (technically) natural have been given. 
This issue is beyond the scope of this paper, but we refer to [29] for a discussion. Secondly 
we have not discussed the thin-shell mechanism so far, which is another essential ingredient 
for the viability of the chameleon mechanism in that it suppresses fifth force modifications 
to e.g. planetary orbits. We will do so in section [6] where radial solutions are investigated in 
some detail. 

3 The disformal/bimetric perspective 

The chameleon mechanism presented in the previous section hinges on the existence of two 
conformally related metrics g^i^ and g^^^, respectively used to construct the gravitational part 
of the action (the Ricci scalar) and the matter part of the action (i.e. matter fields \l/j are 
minimally coupled to g^u). The question we want to answer in this paper is: What happens 
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Figure 1: Left: Plot showing the effective potential Vf.ff{(j)) (solid line), the conformal factor A{(j)) = e^^'^ 



^ (dashed line) in arbitrary units, cf equation (2.6 1. 
= 1 — 6 from bottom to top. Note how the position of 



(dotted line) and a runaway potential V{(j)) ~ (f) 

Right: Plot showing the effective potential for fci 

the minimum (pmin is shifted to smaller (j) and the curvature and hence mass of oscillations around ipmin is 

enhanced as ki increases. 



when this relationship is modified? To be more specific, are there other classes of scalar- 
tensor theories that offer qualitatively distinct implementations of the chameleon mechanism 
(which we take to be an environmentally dependent generation of mass for an otherwise light 
cosmological scalar 0)? 



3.1 A bimetric scalar-tensor theory 

In this section we want to lay out the problem in more general terms without imposing as 
restrictive a metric relation as in section [2j Consequently let us start with a schematic action 
of the following form 

S= d^x^—R + S^ {g^,, ^,) + S^ (3.1) 

where \l/j are matter fields minimally coupled to ^^j, as before and iS,^ denotes an action giving 
the scalar field dynamics of its own. We emphasize that there is no a priori requirement 
constraining S^ to be formed with either g^^y or g^j^^. To enable comparison with the existing 
chameleon literature we fix the form of S^ such that 



S 



d Xy/g [ 



fM^' 



R + X- V{4>) + S^ {g^,, ^, 



(3.2) 



i.e. S^ equips with a canonical kinetic term and a potential minimally coupled to (^^j,. In 



order to investigate (3.2) it now becomes necessary to specify how g^i, and g^^ are related, 
and in particular how enters this relation. For this it will be useful to schematically 



write (3.2) as 



S = I d Xy/gCE {g^,u, 4>) + 



/ d^'x^fgLrn {g,iu, 



^,;1 



(3.3) 



3.2 The metric relation 

In section |2] gravity and matter metrics were conformally related by 

~g^, = A\<P)g^,. (3.4) 

Here we are interested in investigating more general metric relations. A systematic way of 
constructing a more general matter metric in the presence of a gravitational scalar is to 
write 

9,,u = 9,,u{9,,u,(l),d(j),dd(f),...), (3.5) 

which is therefore a function of the gravity metric g^,y, the field (j) and its derivativesF] We 
will truncate this expansion in derivatives at first order, i.e. D^ contributions and other 



higher order terms are ignored. If (3.2) is an effective field theory then subsequently higher 
derivatives of the field can be suppressed by higher powers of the cutoff scale, motivating 
such a truncation. Note, however, that such an argument naturally does not discriminate 
between e.g. (90)" and d"'(j). So in addition, and more heuristically perhaps, we will allow 
operators such as (90)", while disregarding d^(j) for n = 2 in particular. This setup trivially 
protects us from having equations of motion with higher than second order derivatives of 
the field, where we would expect ghost-like instabilities via Ostrogradski's theorem [31] . 

Arguing in a similar fashion, Bekenstein postulated [22] that the most general gravita- 
tional and matter metrics satisfying these conditionq^ give rise to a bimetric theory where 
Qfj^u and Qfj^u are "disformally" related by 

~g^, = A\<f),X)g,, + B\<j),X)d,<l)d,<l), (3.6) 

where we again recognize X as 0's kinetic term. Note that we have implicitly required the 
metric relation to only be a function of coordinate invariants (hence the dependence of A^ 
and B^ on X and not (9^0 by itself). This generalized metric relation therefore consists 
of a conformal piece, with A^ (0, X) being the conformal factor, and a so-called disformal 
piece, where we shall call B^ (0, X) the disformal factor. Interestingly the disformal piece 
of the metric relation exactly mimics the structure of induced metrics in models that are 
intrinsically higher-dimensional in nature. For example, in higher-dimensional dark energy 
setups such as DGP [8] or DBI galileons p3j the induced metric on our 4-brane is g^u = 
Qnu + d^(j)di,(j), i.e. exactly of the disformal type with a trivial disformal factor B^ = 1. 



^While this paper was being finished an analogous construction for a metric that is a composite of several 
other fields (j)i was investigated in [30] , 

^In fact Bekenstein assumes what is essentially a 4D effective field theory in an overall Finslerian geometry 
with gravity and matter metrics related by a single degree of freedom up to first order in its derivatives. 



He then proceeds to show that this has to reduce to a Riemannian geometry described by (3.6). For details 
see [32] , 
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3.3 Generalized equations of motion 



With a matter stress-energy tensor defined as in (2.2) and metrics related by (3.6), we can 



now write down the equations of motion for the general action (3.2). Varying (3.2) we find 




9r9.uu99^^u „ OCe , Idg^^ / 9rf.uA , ^ 9 ri^av^ dg^u 




T^"^ = V,^ + -^V, ( ^/^T^'^l +-^/^r^''V,^. (3.7) 

Note that these expressions neither assume anything about the form of T^'^ nor rely on 
a conformal relationship between g^j^y and g^i, as considered above. It is, however, worth 
emphasizing that T^^^ is explicitly a stress energy tensor not including any contributions 



from S(f, in (|3.1|), so it does not satisfy the Einstein equations by itself. Instead, after having 

T(matter) | rji{(f>) 



mapped all quantities into the Einstein frame, one finds Gf^,^ = SttG IT^^^ '^^ + T^ 



4 Conformal chameleons II: Derivative setups 

Here we wish to focus on the effect of reducing the bimetric relationship (|3.6[) to the simple 



case B{(f),X) = 0. This is the most general purely conformal relation permitted by (3.6). 
The conformal factor ^4^(0, X) is an arbitrary function of the field (p and the first higher 
order coordinate invariant X, (b^s kinetic term, and we have 



~g,, = A'i^,X)g^,. (4.1) 

As such this means the full action under consideration is that of a generalized conformal 
chameleon 

S = j d'x^ (^R + X- V{(P)\ + Sm {A\4>, X)g,y, ^,,) . (4.2) 

The relation between matter stress-energy tensors in different frames straightforwardly gen- 
eralizes to T'^'^ = A~^{(l), X)T^^^ . We will now first compute the associated equation of motion 
for and then comment on the impact a generalized conformal factor has on the chameleon 
mechanism. 

4.1 Equation of motion and effective potential 



The equation of motion for that follows from (4.2) can be written in the following form 



using (3.7) 



-V^ + A^^A^f = -Uct> + Y,J^. (4-3) 



where □ = V^V^^ and 



Ji = -A^xA^'TUcj) 
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J3 = A^f{AA,xx + SA%)U 

^4 = -A=^A,X^;..9>V„^^^ (4.4) 

and we have defined 11 = V^Vj^^V^^V^^ = V ^j,du4>d'^4'd^(t). We have assumed that T^" is 
symmetric and expressed derivatives acting on T^^^ in terms of matter frame variables to 
avoid mixing operators and variables defined for different frames. More explicitly this means 
using 

V^T'^" = V^T'^" - 2rJf^T")^, (4.5) 

where F^^ is the connection associated with transformations between g^y (~ Jordan) and g^y 
(Einstein) frames (for details please see the appendix). Additionally there are two further 



types of terms one might naively expect to arise from varying (4.2) 



Js oc A'A%f^''Ii,pdp<pdP<P 

Je oc A^A^xA^^dpct>d,ct>ff\ (4.6) 

where H^p = V^Vp0 = V ydpcf). However, the symmetry imposed on the conformal factor 
(i.e. A only being a function of coordinate invariants) means that contributions proportional 
to jTs and Jq cancel and hence do not appear in the equation of motion. As an immediate 
consequence there is no direct, Vainshtein-like coupling between the stress-energy tensor 
and derivatives of the field, but matter only enters via TJj Taking g^i, inside its covariant 
derivative and using that Vs = Vs = ds for any scalar s (i.e. covariant derivatives related 
to both metrics act on scalars in the same way), the overall equation of motion can therefore 
be written 

-V^ + A^^A'f= - {l + A^fA^x)^(t) 

+ 2A^f{AA^x4> + ^A^xA^^)X 

+ A^T (AA XX + 3^2^) n 

- A^A^xd''(t)dJ'. (4.7) 

We now proceed to simplify these expressions by considering a uniform matter source, 
making the following key assumption about the system under consideration 

• The stress-energy tensor describes a pressureless, non-relativistic fiuid (Tq = —p) with 
all other stress-energy tensor components vanishing. 

''Note that, at any rate, the Vainshtein-Hke coupUng referred to here would be a derivative coupUng which 
leads to a screening, i.e. a classical renormalization of the kinetic energy, that depends locally on T'^" . This 
means that there will be no screening effect even a small distance away from the source. This is in contrast 
to Vainshtein screening sourced by derivative self-interactions of the scalar 0, e.g. a Xn0 term in the action, 
which will lead to screening inside a Vainshtein radius ry that can extend beyond the source itself, subject 
to introducing an appropriate coupling between matter and <f) such as the linear <f>T. 



Later on (in section pj) we will also assume a static, uniform source in the matter frame 
(Vq.T^'^ = inside the source). Combining this with the first assumption, this is equivalent 
to assuming V^p = dap = (again, inside the source). This means J74 = everywhere 
except in the transition between source and surroundings. For this reason we will keep the 
full J74 term when evaluating radial profiles across this boundary and not make any further 
assumptions about J^ for the time being. Note that we also do not assume a static profile 
for {OqcJ) = 0), which one may want to impose for further simplification in a late-universe 
system like the solar system, which has had time to settle. Modeling matter as a pressureless, 
non-relativistic fluid results in 

Ji = A,xpO(j) 

Js = -p{A^xx + SA-'A%)U 

J4 = A^A xS^S^A-^^p) , (4.8) 

where we have substituted A^T = — p (or equivalently p = A'^p), since p is a conserved 
quantity in the Einstein frame. Explicitly working out J'4, which will be relevant when 
modeling the transition across matter boundaries, one finds 

J4 = 6A-'A^xA,^pX + 3A-'A^^pU + A^xd''(pdaP. (4.9) 

This considerably simplifies the equations of motion, revealing additional symmetries that 
arise due to the functional form of A in combination with requiring a non-relativistic, pres- 
sureless fluid as the matter source. Consequently the equation of motion for (p may be written 
as 

V^ + A^^p = (l-pAx)n0 + 2M,x0^ + M,xxn-Ax5>9,p. (4.10) 

An instructive way to think of the physical properties of this system is to explicitly write 
it as a Klein-Gordon equation with an effective potential that only depends on and "friction 
terms" that encode the dependence on higher derivatives of (j). If we are purely interested in 
the profile inside the source (and hence ignore dap), this means we can write 

Veff,M = □0 + j-i(0,x,p)x + j-2(0,x,p)n 

= n0 + "friction terms" . (4.11) 

To make the form of the effective potential explicit, one can Taylor-expand the conformal 
factor A in powers of X, writing 

A(0,X) = A(°)(0) + A«(0)X + OiX^), (4.12) 

which allows us to write down the effective potential for (p as 

Veff,M = (K0 + M!J^) (1 - A^'^pT' ■ (4.13) 



4.2 Phenomenology and comments 

Given some original potential of the runaway form V{(j)), the chameleon mechanism generates 
an environmentally dependent effective potential that gives a large mass in high density 
regions. Here "high density" and "large mass" are essentially references to solar system con- 
straints [13] on the presence of a fifth force mediated by a scalar degree of freedom. As such, 
any theory with a runaway V{(f)) that successfully implements the chameleon mechanism at 
the very least has to tick two boxes. Firstly it needs to give rise to an environmentally- 
dependent effective potential which has a minimum. And secondly the mass of small oscil- 
lations around that minimum has to be large enough to satisfy fifth force constraints. With 



this in mind let us investigate the effective potential described by (4.10) and (4.13). Note 
that in section |5] we will give explicit examples illustrating each of the effects outlined here 
in detail. 

Position of minimum: For a finite density p and conformal factor A (more precisely, 
A^^^), (4.13) shows that a potential minimum requires 

V;^(0min)+M!?(0m^n)=O, (4.14) 

which is identical to the condition for the minimal conformal chameleon discussed in section[2l 
In other words, the position of the minimum, 0mm, is not altered by the introduction of an 
X-dependent conformal factor. This also shows that a conformal factor which does not 
depend on itself, but only on derivatives of 0, (i.e. A^J = 0, as is the case for A = A{X)) 
cannot generate an effective potential with a minimum, if the original V{(j)) does not already 
possess a minimum itself. The 0-dependence of A is consequently essential to obtaining a 
successful implementation of the chameleon mechanism, if starting with a runaway potential 
(e.g. V{(f)) ~ 0^"). A potential V{(j)) with very small mass my can however be uplifted by 
a pure derivative conformal factor (e.g. A{X)), leading to an effective potential with mass 
m^ ^ my. Below and in the next section we will give explicit examples illustrating this 
behavior. 

Effective mass: A derivative-dependent conformal factor affects the curvature of the 
effective potential. This means derivative chameleons generically come equipped with a 
new mass-altering mechanism. The effective potential, and hence the mass of the field, is 
classically renormalized by introducing higher order invariants into the conformal relation 
as follows 

"^^ = Veff,4><t,{(pmin) = f K0<?i(0™n) + M,°^(0mm) j (l " A'^^\(f)min) p) 

(1-A«(0™„)p)"\ (4.15) 



''standard 



where fn^^^^^^^^ denotes the effective mass for small oscillations around the minimum for a 
theory with identical conformal factor in the limit ^1(0, X — )• 0). In other words, 'ml^^^^^^^ 
is the effective mass for the theory in the limit where higher- derivative contributions can be 
neglected. This mass-altering mechanism can be separated into three branches, which we 
will discuss now. 
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Mass-lifting, ghost-like instabilities and anti-chameleons: Equation (4.15) shows 
that the effective mass m^ is enlarged when < A^^\(l)rnin)p < 1- This is interesting since 
it suggests that derivative chameleon models provide an additional mass-lifting mechanism, 
potentially alleviating the fine-tuning involved in obtaining a sufficiently large mass in dense 
environments for standard chameleon models. However, care must be taken when considering 
mass-lifting solutions for the following reason. Suppose we consider a conformal factor A 
such that a mass-lifting mechanism is in place, i.e. < A^^' ((pmin) Pi < 1 for some given 
energy density pi. Now, assuming A has no density-dependence itself, one can solve for 
a (larger) critical density pcrit above which the solution becomes unstable. The effective 
potential switches sign since 1 — A^^^cpmin) becomes negative, so that (pmin becomes ^max 
and K//,00(0mm) turns negative. Thus we are left with a negative "mass term", signaling 
instabilities, and the solution becomes ghost-like. Figure [2] illustrates these different regimes. 
That different energy densities p will interpolate between stable mass-lifting and ghost-like 
solutions can also be seen from the relevant part of the equation of motion 

V^ + A^^p = (l-M.x) □</> + ..., (4.16) 

where, if pAx > 1, this can be traced back to an action with the "wrong" sign for 0's kinetic 
termo 

There appear to be two obvious solutions to this instability problemjj Firstly one could 
consider making A{(j),X) a function of the energy density p as well. However, especially 
given that the matter stress-energy tensor T^'^ is a variation of the matter Lagrangian with 
respect to the metric g^^ = A'^{(j),X)g^i^, it is not clear what such an iterative dependence 
on T^'^ would mean. Nevertheless it will be an interesting task for the future to think about 
whether there is some convincing way of implementing such a dependence. In any case, 
from a purely phenomenological point of view, a density dependent A allows us to have a 
stable, derivative dependent mass-lifting mechanism for all p. Secondly one may choose a 
conformal factor such that pA^x^^ up to the density cutoff of the theory. Note, however, 
that the effective chameleon mass m"^ is proportional to {pcrit — p)-, where we have defined 
a critical density pcrit such that A^^\(j)min) Pcrit = 1- Introducing a cutoff at pcrit will then 
render derivative-dependent effects on m^ suppressed by that same cutoff scale 

^2 _ ^standard _ 2 f 1 P ] (A 17) 

'"'"/- "N zldUj, N ~ '"'standard \ ^ ^ I ■ \^-^' ) 

[Pcutoff - P>A(^){(pmin> V PcutoffJ 

If the cutoff is low enough, significant derivative-dependent effects on m'^ may still be ob- 
tained, but for a high cutoff density such as the Planck density pp they will be strongly 



^Jumping ahead slightly, also note that, if A,x is a function of (j) and the field is not approximately 
constant (j) ~ (tymin inside the source (the so-called "thick-shell regime - see section [6|, one needs to be aware 
that pA^xi4') being smaller than unity for some initial 0i(r — 0) no longer guarantees that this remains true 
for all values of 4> taken inside the source. 

^Should we decide to bite the bullet and accept the existence of ghost-like solutions above some density- 
scale p, such an approach will also face major challenges when confronted with high energy/density early 
universe physics. 
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suppressed. 

The third branch of mass-altering solutions corresponds to the case when A^^^cpmin) < 0? 
Then the effective mass of is reduced, counter-acting chameleon screening effects 

2 

m" = , "f^ ,,, . (4.18) 

This branch is free of instabilities and provides a robust mechanism to suppress mass terms 
in models with derivative conformal factors A, since the derivative dependence reduces the 
curvature of K//. This suggests that the simplest A^(0, X) models, where A is independent 
of p and no ghost-like instabilities arise for any p, are anti-chameleon models in the sense 
that the effective mass m^ is reduced compared with the non-derivative chameleon limit 
A{(j),X ^0). 

Thin shell regime: The computation of Veff and its minima/maxima and masses above 



was oblivious to "friction" terms in (4.11) (by definition, since the effective potential is a 
function of (p and not its derivatives). But such terms automatically arise as a consequence 
of a higher order conformal coupling. Here "friction terms" is a reference to terms other than 
0(f) that have a derivative dependence on 0, e.g. JJ2 and J'3 which encode the dependence on 
X and n. While not influencing Veff, these terms do impact the dynamics of 0, warranting 
further investigation. 

One particularly interesting consequence is a modification of the so-called thin-shell 
regime. In typical chameleon theories the chameleon-charge of large, massive objects can 
be Yukawa-screened with only a thin-shell on the outside of the object contributing to the 
exterior 0-profile [3111]. This is essential for e.g. avoiding unacceptably large effects on plan- 
etary orbits due to a chameleonic fifth force. Now the presence of additional friction terms 
modifies the gradient of inside the source and one therefore expects a modification of the 
thin-shell effect as well. We will discuss this in detail in section |6| where we investigate radial 
solutions around massive sources, focusing on new-found phenomenology due to derivative 
conformal factors. 

Equivalence principle (violations): A final comment on possible equivalence princi- 
ple violations. By construction, chameleon models (both standard as well as the derivative 
generalization considered here) respect the weak equivalence principle. In the non-derivative 
case, the extra degree of freedom locally influences the dynamics, but it does not discrimi- 
nate between different test masses/types of matter due to the universal coupling of ^4^(0) to 
all matter fields \E'i (field-dependent couplings are discussed in [S])- In the derivative case, 
an additional degree of freedom X enters. But since all matter still universally couples to 
A'^{(j),X)gfj^i,, all test masses locally experience the same gravitational force (note that is 
viewed as a gravitational scalar here). In both cases 0's profile of course does not remain 
constant across space-time, as it depends on the ambient density. This means the strong 
equivalence principle is trivially violated (as is, in fact, the Einstein equivalence principle 
- cf. with e.g. the Horndeski constructions in [3S] where the strong, but not the Einstein 
equivalence principle is broken). 
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The main point we wish to make here is that a prima facie worry one might have when 
considering yl^(0, X), namely that a derivative-dependent couphng to matter violates even 
the weak equivalence principle, is not justified. This is simply due to the fact that there is 
no dependence of the gravitational coupling on the momentum of matter (\E'j) test masses, 
but only a derivative coupling to the gravitational scalar 0. Nevertheless this is an area that 
warrants further investigation, since known violations of the strong equivalence principle 
in chameleon models (cf. [361 EZl EHl EHl SO]) will be modified by introducing a derivative 
dependence. Computing this in detail should enable further disentangling of derivative vs. 
non-derivative chameleons. 

5 Effective potentials and the chameleon mass 

In this section we illustrate how chameleon mechanisms arise in derivative theories with 
a number of explicit examples. We focus on the mass-altering mechanism for different 
conformal factors. 

5.1 Example I: Taylor-expanding A{(j), X) 



We begin by reminding ourselves of (4.12) and Taylor-expanding the conformal factor 



A(0,X) = A(°)(0) + A^'\<P)X + OiX^), (5.1^ 



where the associated mass of oscillations around an effective minimum is given by (4.13). 
We will ignore higher orders in X and, as a first simple example, focus on the zeroth order 
contribution in to A^^\ treating it as a constant. As such the conformal factor takes on 
the form 

A{(j),X) = A^''\<f)) + kxX. (5.2) 

The resulting potential and effective mass can be computed to give 

-L ^ rvxP 

1 - kxp 1 - kxp 

This illustrates the point made in the previous section about the existence of a mass-raising 
and a mass-lowering branch for derivative chameleon models. In the simple setup considered 
here, if kx is negative, the chameleon mass m^ is reduced compared with the standard 
A = A(0) theory. The fine-tuning necessary in order to get a sufficiently large chameleon 
mass is therefore made more severe in this case, since an additional mass-lowering mechanism 
is at work. For positive kx a mass-lifting mechanism operates, increasing m'^. However, for 
kxP = 1 the mass diverges and, once kxP exceeds unity, the "mass" turns negative and 
the solution ghost-like. As discussed in the previous section, this is cause for concern, since 

13 



given some gravitational theory the conformal factor is fully specified (and hence kx is fixed 
in this case). As such, sources with densities above pcrit = l/^x essentially see an unstable 
inverted potential —Vefj inside the source. The bottom left graph of figure [2] illustrates 
this point by plotting the dependence of the effective potential on kx for some given source 
density pi. This amounts to considering different normalizations of the effective potential. 
On the other hand, the bottom right plot shows how varying p for some given, positive and 
fixed kx affects the solution. Note that having restricted to positive kx means there are no 



mass-lowering solutions present in this plot. As expected from (5.3), changing p modifies 
both the position of the minimum as well as as changing the potential's normalization and 
hence curvature/mass via the 1 — k^p term. 

As an aside, notice that, if one does allow the conformal factor to depend on energy 
density p (for a discussion of this approach see the previous section), then one can construct 
a solution which remains ghost-free inside sources with arbitrarily large energy densities p. 
E.g. we may impose 

/l(0, X) = e'^^^ + p-^ (1 - e'^^^) X + 0(X2), (5.5) 

which results in an effective mass given by 

^2 _ gfe2<^_ {V^^{(l)min) + fc^pe'^^*— ) . (5.6) 

In other words, in this particular example /c2 allows tuning the mass m'^ arbitrarily. The 
potential is stable for all and the mass can be altered by modifying the X-dependence of 
the conformal factor A^^ 

5.2 Example II: A separable A{(f), X) 

In the Taylor-expanded picture laid out in the previous section, considering a separable 
conformal factor A(0, X) amounts to setting A*^°)(0) = A^^^cp). An interesting feature 
compared with the kx case discussed above is that the normalization of the potential now 
also becomes a function of 0. To see this explicitly we write the conformal factor as 

^(0,X) = i?(0)C(X), (5.7) 

which means the effective potential satisfies 

Veff,M = {V^ + pB^^C^"^) (1 - BC'^'^pY' , (5.8) 

where C* refers to the i-th component in a Taylor-series expansion of C{X) in powers of X*. 
Since at the minimum of the potential Vef j ^{(pmin) = 0, we obtain an effective chameleon 
mass for oscillations around (pmin of 

"" " 1 - \k2B{ct>m.n)Ci^)p • ^ ^^ 



^If one is concerned that additional singularities are introduced by the p ^ dependence of A here, notice 
that a condition such as A'-^'p ^ 1 — e~^ will also ensure a ghost-free mass-lifting regime. 
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Viewing the effective potential as a renornialized version of tlie corresponding standard 
(i.e. non-derivative) chameleon setup with C = 1, the effect of the B{(j)rnin) term in the 
denominator can be described as making this renormalization 0-dependent. The top right 
plot in figure [2] shows the effective potential for a separable conformal factor of the form 
A{(j),X) = Exp[ki(f) + k2X]. Note that we require a positive k2, because otherwise the 
conformal factor diverges as d^(j) ~^ 0? i-^- iio stable, static solution is possible. Now the 
effective mass can be written as 

^, = K00(0^.n) +Pk!^'''- = ^2 (1 _ fc,e^-^-v)-^ . (5.10) 

The solution again separates into the three branches discussed. For positive -8(0) and when 
< ^2-8/) < 1 the effective potential has a well behaved minimum protected by an infinite 



potential barrier at (l)crit-, where (^crit satisfies k2B{(f)crit)p = 1 for a given p. Equation (5.10) 
then shows that m'^ is enhanced by the derivative-dependence of A{(j),X). For values of 
> <Pcrit there do exist unstable regions of parameter space. However, a particle, which 
starts at an initial field value 0j for which A;2-B(0i)p < 1, is protected from entering the 
unstable region. Conversely, given a particular conformal factor that fixes k2, solutions with 
positive B{(j)) inside a source with energy density such that p > (^2-8)"^ are unstable. Finally, 
for negative B{(j)) a stable, mass-lowering solution is obtained similar to the one discussed 
in the previous section. The top right plot of figure |2] illustrates these three branches and 
one can see the generation of an effective minimum and how the mass m^ can be tuned by 
varying /c2 (for a fixed p) at the expense of lowering the critical field value (pcrit- 

5.3 Example III: A purely derivative conformal factor A{X) 

Suppose we have a purely derivative conformal factor that does not depend on the field value 
of itself. This is interesting, since it means the theory has a shift symmetry —)■</) + c, 
potentially protecting it from a number of quantum correctionaj (assuming this is at most 
softly broken by the potential V{(j)) - cf. [H], where this point is discussed in an inflationary 
setting). As such, let us consider a conformal factor of the following form 

~g^, = A\X)g,, = (1 + A(i)X + ^(X^)) g^,, (5.11) 

where A^^^ has been appropriately normalized (i.e. there is no fundamental reason why 
it should be unity) and where all A^ are now constants and therefore not functions of 0. 
The reason to require A^^' 7^ is that otherwise any coupling between gravity and matter 
vanishes once 9^0 = 0, e.g. once has settled into its minimum. The full action can then 
be written as 

S= I d^x^ (^R + X- 1/(0)") + Srn {A\X)g^,, ^,) . (5.12) 



^This statement is of course not as strong as for a theory with an additional Galilean shift symmetry 
(cf. [9]) where an effective non-renormalization theorem can be derived |4Tl|42]. Note, however, that even 
there analogous concerns enter due to the non shift-symmetric matter coupling (/)T, cf. discussions in [43 [ 144 ) . 
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Figure 2: Top Left: Plot showing the effective potential Vef /{(/)) (solid line) for a separable conformal 



factor A{(I),X) = B(0)C(X) 



— pki4'+k2X 



(dashed line) and a runaway potential V{(f)) 



(dotted line) 



in arbitrary units. Top Right: Logarithmic plot showing Veff{4>) for different choices of fc2 (negative, zero 
and positive from bottom to top with p = 1 in arbitrary units). This shows how the amplitude, and hence 
mass, of the chameleon field grows as fc2 is enlarged, while the stable region of parameter space is reduced. 
The dashed line shows how the pole in the denominator of Veff creates a discontinuity at (jjcrit and the 



effective potential assumes a runaway form for values of (/)>q 



In other words, the solution transitions 



between the mass-lifting and ghost-like branches at 4'crit- Note that the value taken by (l)crit is a function 
of p and the effective potential. Bottom Left: Plot showing the effective potential for a conformal factor 



as in equation (5.2) and for different values of the parameter kx and fixed p: The solid line shows the non- 
derivative chameleon with kx — 0, dashed lines show the mass-lifting branch with < kxp < 1, dot-dashed 
lines show the mass lowering branch with fc^ < and dotted lines lie in the unstable region with kxp > 1. 
Bottom Right: Analogous plot varying p for some fixed positive kx (hence no mass- lowering branch is 
present here). Varying p changes the position of the minimum as well as the normalization of V^ff and 
hence its curvature and mass m^. p is zero for the solid line and below/above Pcrit for dashed/dotted lines 
respectively. 



Now we have already seen that a conformal factor, which does not depend on the field 
value itself, cannot create an effective potential with a minimum, if V{(j)) is of the runaway 
form (so it does not have a mass term). But what if dark energy is in fact sourced by a very 
light, but not massless, cosmological scalar? Then V{(f)) does already have a (very small) 
mass term and hence a minimum. Let us focus on a particularly simple toy model and 
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consider a power-law potential of the form 

V^*) = \{£)\ (5.13) 

i.e. a simple mass term for m^ = 0q^^, which allows us to tune the mass of the field by 
choosing 0o- In order for this field to act as dark energy on large scales, it has to be 
extremely light m^ = 0q ^ ^ -f^o ~ 10^'^^ey. Now the effective equation of motion for is 

V^ = {1- M.x) □</> + pA,xxn - A^xd^cPVaP- (5.14) 

and we can write down an effective potential Ve//(0) = V{(f)) (l — A^^^pj . In dense envi- 
ronments the effective mass of the chameleon field consequently goes as 

^' = Y^^^ = {^i-^iA'''py\ (5.15) 

which can be very large subject to < 1 — A*-^^p <^ 1. This outlines how a mass-lifting mech- 
anism caused by a purely X-dependent conformal factor can give rise to a viable chameleon- 
type solution. An otherwise extremely light cosmological scalar field then acquires a large 
mass in dense environments. 

We emphasize that the point made here is independent of the exact form of the potential. 
Given any potential for a non-massless field 0, i.e. a potential with a minimum at (pc, the 
mechanism outlined here will raise the field's mass in a density-dependent manner. We can 
therefore apply this purely derivative mass-lifting mechanism to any potential V^(0) for a 
field with a small mass m^ ^ 0, which reproduces the desired dark energy behavior on 
cosmological scales. 

However, this approach faces a major obstacle. The field's mass has to be enlarged by 
several orders of magnitude if a light cosmological scalar is to escape detection on solar 
system scales. This means A^^^p ~ 1 for a large range of densities p, so that A^^^ has to depend 
on p in order to suppress the otherwise strong dependence of m} on p and in order to prevent 
ghost-like instabilities from developing. We refer to section |4] and the previous examples 
in this section for a discussion of the possibility of such a density-dependent conformal 
factor. Also note that the mass-lowering branch with A^^^ < 0, albeit perhaps less interesting 
phenomenologically, does not face this problem just as discussed in the previous examples. 

6 Radial solutions, the thin-shell effect and "friction 
terms" 

In this section we investigate the static, radial 0-profile in and around a spherically symmetric 
body with uniform density p - a good approximation for the profile around the sun or earth, 
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for example j^ Expressing the general equation of motion (4.10) for such a profile one finds 



V.^ + A., 



:i-M, 



X, 



(Pep 2 dcp 
dr'^ r dr 



P^,X4> 



dr 



+ pA 



d^ 



XX 



dr"^ \ dr 



- A^xd'^^d^p, 



(6.1: 



where p and (p are functions of r. We now modify the approach taken in |4] in that we still 

assume 

Pc for r ^ i? 

Poo for r <^ R. 
where R is the radius of the spherical object, pc is its density and p^ 



p{r) 



(6.2) 
is the density of the 



surroundings. With an eye on computing the gradient term dap{r) in (6.1), we model the 
boundary between source and surroundings by a very sharp, but smooth, transition between 
Pc and poo- The particular template we shall adopt is 

1 



P(^) = 7^{Pc - Poo)(l - Tanh [s(r - R)]) + p^ 



(6.3) 



taking s ^ 1, so that p effectively remains at its asymptotic values except for a sharp 
transition around r = R. A unique solution for (6.1) requires specifying two boundary 
conditions. Again following |1], we take these to be d(j){r = 0)/dr = 0, so that the solution is 
non-singular at the origin, and — )■ 0oo as r — )■ 00, which ensures that the 0-mediated force 
between two test bodies vanishes as r — )■ 00. 

In what follows we will first recap how this setup gives rise to a thin-shell effect for the 
minimal conformal chameleon (i.e. the standard case). Then we discuss how radial solutions 
and thin-shell behavior are modified by introducing derivative conformal factors, paying 



special attention to the effect of "friction terms" in (6.1) (i.e. X and 11 dependent terms as 



mj2 and J3 (4.8)). 



6.1 The thin-shell effect for standard chameleons 

Any given physical system will come equipped with a specific V^(0), p, R etc. This will 



then allow us to compute the corresponding solution to (6.1), in particular fixing the initial 
field-value 0i at r = 0. Alternatively one can explore the system's solutions by choosing 
and finding the corresponding region in {p, R, ...) parameter space |l]p 



some 



Now for the minimal conformal chameleon, as described in section |2| the radial equation 
of motion is simply 

2: 



Vs + A„ 



/d^ 

\ dr"^ r dr 



(6.4) 



®For a time-dependent chameleon field in the case of a radially pulsating mass, see [46] . 

^This amounts to solving (6.1) as a classical mechanics problem with <j> being a position- and r being a 



time-coordinate. Note that the opposite signs for temporal and spatial dimensions in g^xjy's signature mean 
that, for the radial solution we are considering here, the particle whose "position" is described by (p "moves" 
on the inverted, and hence unstable, potential —Veff- 



^ 



50c 



-# 




Figure 3: Left: The radial solution for a standard chameleon with conformal factor A{(f)) — e*^^"^ and 
different values for the parameter fci = 1,1.05,1.1,1.2 from bottom to top. The initial 0i(r = 0) remains 
fixed for all fci here (note that the position of the extrcmum (pc also depends on A and hence fci). Right: 
Plot showing how a suitably chosen (jji leads to a radial solution that asymptotes to (/)oo as r — >■ oo. fci = 1 
here. 



Outside the source, i.e. at a radius r > R, the solution for then assumes the form jl] 



(p[r) ~ 
AR 



AR\ M,e^™-('^-^) 



R 



R — R 



roll 



R 



$. 



R 



<1, 



(6.5) 



where R is the radius of the source, Mc is its mass, 0oo and moo are the field value and mass 
of the field as r — )■ oo, 0c is the local extremum of the 0-potential inside the sourcq^and $c 
is the Newtonian potential at the surface of the source. We have also assumed that AR/R, 
the so-called thin-shell suppression factor, is small. An intuitive picture for this setup is the 
following: The field is released at "time" r = from 0j. Inside the source the field remains 
stuck near (pi with its dynamics dominated by the 2/r ■ d(f)/dr friction term. Eventually, at 
r = Rroih the field starts rolling and hence developing gradient terms. Only contributions 
from the region between Rroii and i?, i.e. the region where gradient terms develop, are felt 
by the exterior profile. Since AR/R <^ 1, the chameleon force from a large massive body 
on a test mass is therefore thin-shell suppressed. This is essential in explaining why e.g. 
planetary orbits are not affected by a chameleonic force. In the language set out above the 
thin-shell regime corresponds to 0j — 0c ^ 0c- 

In contrast, the radial solution for the unsuppressed "thick-shell" regime, where 0j>0c, 
is given by 



Mg-m«>(r-/?) 



r 



(6.6) 



^"^Since we are effectively dealing with motion along —14//, (f>c is the same as (jimin in previous sections. 
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Here the field starts significantly displaced from 0^ and consequently starts rolling almost 
immediately. Gradient terms develop inside the source and there is no thin-shell suppression 
(hence the absence of the AR/R factor). It is intuitively clear that the presence of further 



derivative terms in the general solution (4.10) will modify this behavior, since it will affect 
gradient terms in (p. We now move on to explore such "friction" terms and their effect on 
the thin-shell mechanismP^ 



6.2 Radial solutions for derivative chameleons 



5*1 






Figure 4: Left: Plot of the radial solution for a conformal factor as in (5.2) with kx — (solid line), kx 



increasingly positive (dashed lines - kx increases as one moves out from the solid line) and kx increasingly 
negative (dotted lines - kx decreases as one moves out from the solid line). Middle: Plot of the radial 
solution for a conformal factor as in (6.9). Again the solid line denotes kx = 0. Dashed lines are for 



fixed positive kx and outwardly increasing ^2, while dotted lines are for negative kx, again with outwardly 
increasing fc2. In order to illustrate the effect of derivative contributions more clearly the solution is shown 
for an initial displacement (j)i — (pc ^ O((/)c/10), i.e. an intermediate region between thin- and thick-shell 
screening. Notice that the most significant effect comes from the dp term at r = _R, which produces a "kink" 



in the radial solution (actually smooth due to the density profile (6.3)). Right: Analogous plot with the 
solid line describing kx — 0. kx is positive above and negative below the solid line. Here dashed lines 



represent the solution including the friction term that depends on A,xct> (cf equation (6.10)), while dotted 
lines ignore this contribution. The additional "friction term" therefore drives the solution further away from 
the non-derivative limit. 



Derivative solutions without "friction terms": How does the radial solution and 
thin-shell suppression change once we go beyond the simplest chameleon setups? We begin 
by considering a simple conformal factor of the form 

A{ct>,X) = e^^^ + KX, (6.7) 

i.e. a sum-separable A{(l),X) = B{(f)) + C{X) with purely linear dependence on X. This 



reduces (6.1) to 



V^ + A„ 



'I- -A ) t^^l^ 
dr"^ r dr 



- A,x5>a„p. 



(6.8) 



^^ Further discussions of standard thin-shell screening and spherically symmetric solutions can be found 
e.g. in [amillllllS]. 
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In contrast with the standard chameleon case there are therefore two new effects. Firstly 
the (1 — pA^x) term, which effectively renormalizes the potential and is responsible for the 
mass-changing mechanism discussed in section |5j And secondly the new gradient term in 
daP, which will only be significant very close to r = R. Importantly we have no X- or 
Il-dependent friction terms for a conformal factor of the form considered here. 

We plot radial solutions for (6.8) in the left panel of figure |4] for different values of k^- 
The behavior observed can be understood as follows. For positive kx, the (1 — pA^x) term 
increases the curvature of the effective potential (and hence the mass of oscillations around 
the minimum). This means that the driving term K//,0 is enhanced and will overcome the 
2/r-d(f)/dr friction "earlier" in the evolution, i.e. it will reduce Rroii- In other words, in order 
to obtain the same thin-shell screened exterior solution one now needs to release the particle 
(j) from even closer to the minimum value (pc, so that (p'^^'^^'""-^^'"^ <^ ^standard ^ The dap term in 
fact further increases this tendency, giving an additional positive "kick" to the gradient of (f). 
Note that this second effect is oblivious to the profile inside the source and as such does not 
modify the thin-shell condition. However, in order to reach the same boundary value 0oo as 
r — 7- oo, both new terms require an initial (pi closer to (p^ than in the non-derivative case. 

For negative kx the converse is true. The curvature and associated mass of the effective 
potential is reduced and thin-shell screening is enhanced, i.e. Rrou is pushed closer to R. 
This broadens the thin-shell screened parameter-space and reaching the boundary value 0oo 
as r — )■ oo requires an initial 0j further away from 0c than in the non-derivative case for 
negative kx- Overall the new derivative dependent effects found here modify the parameter- 
space for thin-shell screened solutions, with a positive/negative kx weakening/strengthening 
thin-shell screening respectively. As intuitively expected from modifying the curvature of the 
effective potential, the mass-lifting branch is therefore associated with a suppressed thin-shell 
mechanism, whereas the mass-lowering branch enhances thin-shell screening. 

Derivative solutions with "friction terms": Here we finally wish to examine the 
effect of "friction terms" that depend on derivatives of 0. The toy model we adopt has a 
conformal factor 

/l(0, X) = e'^i'^ + kxe^^'^X + C(X2), (6.9) 

where we will ignore the higher order corrections 0{X'^). This means the radial solution ( |6.1 ) 
simplifies as H-dependent contributions drop out (such contributions only come in for a 
conformal order with non-linear dependence on X). The equation of motion consequently 
becomes 

V, + Ayp = (1 - pA^x) (^ + If^ - pA.x, {^ - Axd^'PdaP. (6.10) 

In the middle and right panel of figure |4] we plot solutions for constant ki, but varying ^2 and 
kxi essentially tuning the contribution of the 2pA^x4>X term. One can see that a large /c2 in 
combination with the presence of the daP term magnifies the "kick" at r ~ i?. Perhaps more 
interestingly, for positive kx the direct coupling of to its derivative in the conformal factor 
again counteracts the thin-shell effect, by driving the field away from near its minimum 0c. 
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As before this results in a reduced Rroii- Also as before the converse is true for negative kx- 
The inclusion of the 2pA^X(i)X term therefore strengthens the tendencies observed above: 
The mass-lifting branch becomes yet more fine-tuned in order to obtain a thin-shell screened 
solution, while for the mass-lowering branch the parameter-space corresponding to thin-shell 
screening is broadened. 

This concludes our brief survey of how radial solutions around a massive source are 
modified by the introduction of derivative terms into the conformal factor. As a generic 
feature, we observe that a positive kx, linear dependence on X, which leads to a mass-lifting 
mechanism, simultaneously tightens thin-shell screening constraints, requiring a value of (pi 
closer to 0c in order to maintain a nearly constant field value inside the source. The converse 
is true for negative kx- 

7 Disformal chameleons: A no-go theorem 

Until here we have only been considering conformal couplings to matter. Here we wish 



to investigate whether switching on the disformal B^ term (3.6) can result in interesting 
modifications to the chameleon mechanism (for other disformal dark energy models we refer 
to lini [50]). As such, matter now couples to the full disformal matter metric 

g^, = A\<f),X)g^, + B\<P,X)d^<pdA- (7.1) 

Note that such a metric always has a well-defined inverse [32] given by 

r^" = A-2(0, X) [g^, - B\(t>, X)C-\<I>, X)dyd''(l)) , (7.2) 

subject to the condition that A\(j),X) ^0^ C7(0,X), where C = A\(j),X) + B^{(j),X)X. 
We then have 

r //Vf2 1 \ 

S= I d^x^ ( —R - 2^ - Vi(l)) j + Sra (^2(0, X)g^, + B\<P, X)5^09,0, ^,) . (7.3) 

Just as in the conformal case we consider a non-relativistic, pressureless source. In addition 
we here focus on static solutions for which doip = (requiring spherical symmetry could 
reduce this further to = 0(r)). As such our assumption list now is 

• The stress-energy tensor describes a pressureless, non-relativistic fluid (Tq = —p) with 
all other stress-energy tensor components vanishing. 

• 0's dynamics are described by a static solution (9o0 = 0). 

As before we may also assume a static, uniform source in the matter frame (VqT'^'^ = 
inside the source), but this won't be necessary for the argument here. Instead we are solely 
interested in whether a chameleon mechanism can arise inside a source as a consequence of 
disformal contributions here at all. 
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The equation of motion (3.7) then immediately gives rise to a no-go theorem. This 
is because all disformal (S^-dependent) terms in fact involve a contraction of one of the 
following forms 

a^05,0T''^ 5^0T^^ a^a^^T'^^ d^ctN^f^^". (7.4) 

As such all disformal contributions vanish for a static solution around a non-relativistic, 
pressureless source, since the only non-zero component of the stress energy tensor is T^*^, 
which vanishes when contracted with OqcJ) for a static solution. Chameleonic effects from 
the disformal coupling are therefore suppressed when considering e.g. solutions around the 
earth. We do expect disformal effects to play an important role when computing relativistic 
or non-static corrections to this solution, however. Nevertheless, this shows that, with the 
mild assumptions implemented above, the dominant contribution to the chameleon profile 
necessarily has to come from a conformally coupled matter metric. Disformal effects here 
only come in at next order in relativistic corrections and for non-static profiles. Note that 
we have focused on chameleonic behavior here. As comparison with the discussion of jTs and 
Jq (4.6) in section H^ shows, Vainshtein-like screening can arise as a consequence of disformal 



couplings 
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8 Conclusions 

The main results of this paper can be summarized as follows. 



A derivative chameleon model described by (4.2), where all matter is universally and 
minimally coupled to a metric g^^i, = A'^{(f),X), generically comes with a new mass- 
altering mechanism. This separates into three branches depending on the local energy 
density p and the conformal factor A: A mass-lifting ("screening") branch, where 
derivative effects lead to an additional enhancement of the effective chameleon mass 
compared with the non- derivative A'^{(f), X — )• 0) limit. For high energy-densities above 
some pcrit this transitions into a ghost-like branch in which the effective potential 
becomes unstable. Thirdly there is a stable mass-lowering branch. 

This suggests that derivative chameleon models which are not plagued by ghost-like 
instabilities typically lower the effective chameleon mass. Implementing a mass-lifting 
mechanism while avoiding ghost-like instabilities for a range of different densities p 
requires some additional "engineering" , e.g. introducing either an appropriate energy 
density cutoff or a density-dependent conformal factor. 



^^Once the metric g^,y is promoted to a building block for fm-ther gravitational parts of the action as in [33] 
- e.g. a cosmological-constant-like term ^/gA., an extrinsic curvature K^i, and a Ricci scalar R for g^t^ - then 
a full galileon/Horndeski type solution can be obtained, making the Vainshtein screening in question more 
robust. 

23 



• A very light, but massive, cosmo logical scalar can in principle be chameleon-screened 
by a purely derivative-dependent conformal factor A^(X) in the mass-lifting branch, 
which offers added protection from quantum corrections over A^(0) due to the presence 
of the shift symmetry — )■ + c. 

• The position of the minimum (pmin of an effective chameleon potential K// cannot be 
affected by derivative dependent terms (X-dependent as well as for higher order terms). 
If we start with a runaway potential, e.g. V{(j)) ~ 0"" for some positive n, this means a 
0-dependent conformal factor is necessary in order to produce chameleon screening. In 
other words, in this particular setup A'^{(f),X) can yield chameleon screening, whereas 
A'^{X) cannot. 

• Radial solutions and the thin-shell mechanism are modified in derivative chameleon 
models. The region of parameter-space exhibiting thin-shell screening is reduced in 
the mass-lifting branch and enhanced in the mass-lowering branch of solutions. 

• Disformal contributions to the matter metric cannot source chameleon phenomenology 
for a static solution around a non-relativistic, pressureless source (they can source 
Vainshtein-screened solutions though). 

To conclude, chameleon models provide a well-established framework for reconciling the 
presence of a light cosmological scalar that drives cosmic acceleration with small scale fifth 
force constraints. Here we have shown how chameleon-type setups can be generalized to 
derivative-dependent conformal couplings, suggesting that this naturally generates a further 
mass-changing mechanism. The mass-lifting branch offers the possibility to alleviate the 
fine-tuning involved in making a chameleon fit fifth force constraints and also comes with 
the exciting promise of making chameleon models more robust from a quantum perspective 
by endowing them with a shift symmetry. However, this branch also faces a ghost-problem 
which may be answered by introducing a p-dependent conformal factor or a p-dependent 
cutoff for the theory. For the mass-lowering branch no such problems exist, suggesting that 
a mass-suppressing mechanism is in fact a typical feature of derivative chameleon models. 
As such we hope that this work contributes to the enterprise of providing a wider survey of 
the ways in which chameleon-phenomenology might be implemented, realizing what types of 
challenges it has to face, and therefore putting such models on a firmer footing by showing 
what requirements stable chameleon models have to meet. 
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A Appendix: Covariant derivatives 

Here we present some details of the calculation for a generalized equation of motion in 
conformal derivative chameleon setups, expanding on the results presented in section |4} We 
start with a metric relation 

~g^, = A\(t),X)g,,. (A.l) 

The covariant derivatives for the two metrics are related by 

Vc(^p = Va(^l3 - ^l/s^^y (A. 2) 

For the Einstein frame covariant derivative Vq acting on the matter frame stress-energy 
tensor, one therefore finds 

Vaf^" = Vaf^" - T^paT^" - r^^T'^^. (A.3) 

In calculating the connection F^^ explicitly, the following relations will be useful 

VaA = A^^da(f) + A^xdc,X, 

= A^^da<p- A^xdad^<pd^'<p. 

WaA^x = A^x<i>da(p + A^xxdcX, 

= A,xc^d^<f) - A,xxdc.d^(j)dy. (A.4) 

As such we can work out the connection, arriving at 

= A~^{25l^V^)A-g^pg~<'VsA) 

= 2A^^A-Hl^d^^<p - 2A^xA-'5lpd^)d^<i)d^<i) 

- gapg^'ds<pA-'A^^ + g^pg^'dsd^(PdyA-'A,x. (A.5) 
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